      SUBROUTINE SETUP
      RETURN
      END
C   HS72
      SUBROUTINE EF(X,FX)
      INCLUDE 'FUCO.INC'
      DOUBLE PRECISION X(NX),FX
      ICF=ICF+1
      FX=1.D0+X(1)+X(2)+X(3)+X(4)
      RETURN
      END
      SUBROUTINE EGRADF(X,GRADF)
      INCLUDE 'FUCO.INC'
      INTEGER J
      DOUBLE PRECISION X(NX),GRADF(NX)
      ICGF=ICGF+1
      DO      100      J=1,4
      GRADF(J)= 1.D0
 100  CONTINUE
      RETURN
      END
      SUBROUTINE EH(I,X,HXI)
      INCLUDE 'FUCO.INC'
      INTEGER I
      DOUBLE PRECISION X(NX),HXI
      CRES(I)=CRES(I)+1
      RETURN
      END
      SUBROUTINE EGRADH(I,X,GRADHI)
      INCLUDE 'FUCO.INC'
      INTEGER I
      DOUBLE PRECISION X(NX),GRADHI(NX)
      CGRES(I)=CGRES(I)+1
      RETURN
      END
      SUBROUTINE EG(I,X,GXI)
      INCLUDE 'FUCO.INC'
      INTEGER I 
      DOUBLE PRECISION X(NX),GXI
      IF ( GUNIT(1,I+NH) .EQ. -1 )CRES(I+NH)=CRES(I+NH)+1      
      IF(I .GT. 2)      GOTO 100
      IF(I .EQ. 2)      GOTO  50
      GXI=.0401D0-4.D0/X(1)-2.25D0/X(2)-1.D0/X(3)-.25D0/X(4)
      RETURN
  50  CONTINUE
      GXI=.010085D0-.16D0/X(1)-.36D0/X(2)-.64D0/X(3)-.64D0/X(4)
      RETURN
  100 CONTINUE
      IF(I .GT. 6) GXI=((11.D0-DBLE(I))*1.D5-X(I-6))
      IF(I .LE. 6) GXI=(X(I- 2)-1.D-3)
      RETURN
      END
      SUBROUTINE EGRADG(I,X,GRADGI)
      INCLUDE 'FUCO.INC'
      INTEGER I,J
      DOUBLE PRECISION X(NX),GRADGI(NX),C(4,2)
      SAVE C
      DATA C/4.D0,2.25D0,1.D0,.25D0,.16D0,.36D0,.64D0,.64D0/
      IF ( GUNIT(1,I+NH) .EQ. 1 )RETURN
      CGRES(I+NH)=CGRES(I+NH)+1
      IF(I .GT. 2)    GOTO 200
      DO     100      J=1,4
      GRADGI(J)=C(J,I)/X(J)**2
  100 CONTINUE
      RETURN
  200 CONTINUE
      DO      300      J=1,4
      GRADGI(J)=0.D0
  300 CONTINUE
      IF(I .GT. 6)        GRADGI(I-6)=-1.D0
      IF(I .LE. 6)        GRADGI(I-2)=1.D0
      RETURN
      END
      BLOCK DATA
      INCLUDE 'BLOCO.INC'
      INTEGER I,J
      DATA NAME/'HS72'/
      DATA (X(I),I=1,4)/4*1.D0/
      DATA N/ 4/ , NH/0/ , NG/10/
      DATA DEL0/1.00D0/ ,TAU0/1.0D0/ ,TAU/.1D0/
      DATA ((GUNIT(I,J),I=1,3),J=0,10)/-1,0,0,-1,0,0,-1,0,0,
     F           1,1,1,1,2,1,1,3,1,1,4,1,
     F           1,1,-1,1,2,-1,1,3,-1,1,4,-1/
      END
